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Abstract 

In this paper, we build a Two-Scale Macro-Micro decomposition of the Vlasov equation with a 
strong magnetic field. This consists in writing the solution of this equation as a sum of two oscillating 
functions with circonscribed oscillations. The first of these functions has a shape which is close to 
the shape of the Two-Scale limit of the solution and the second one is a correction built to offset this 
imposed shape. The aim of such a decomposition is to be the starting point for the construction of 
Two-Scale Asymptotic-Preserving Schemes. 

1 Introduction 

The goal of this paper is to make the first step towards the setting out of a new class of numerical 
methods: the Two-Scale Asymptotic-Preserving Schemes or TSAPS. We intend to use these methods 
in order to treat problems involving the following two numerical difficulties: first, dealing efficiently on 
long time scales with solutions having high frequency oscillations (the Two-Scale approach) and second, 
an accurate and stable managing of the transition between different regimes, using a unified model (the 
Asymptotic-Preserving approach). These problems naturally occur in solving a Vlasov equation implying 
some small parameter, for the modelling of the dynamics of charged particles in the presence of a strong 
magnetic field. More precisely, we are interested to develop a model whose discretization will be able to 
simulate both, the regime when the parameter is not small (i.e. the magnetic field is not large), and the 
limit regime obtained when the parameter is small. The discrete scheme will automatically shift from one 
regime to the other. In addition, in the limit of the small parameter, the particles mean behaviour will be 
efficiently described by the two-scale limit while the micro behaviour will be expressed by some correctors. 

To make the purpose more precise, we first refer to [4, 8, 13, 14, 15, 16] and the references therein 
as previous works about different regimes, by taking limits of Vlasov or Vlasov-Poisson equations cor- 
responding to small parameters, having in common the framework of a large magnetic field. Then, we 
refer to [1, 2, 9, 11, 21, 22, 23] for a detailed description of Two-Scale Numerical Methods. As for the 
Asymptotic-Preserving schemes, although there is an abundant literature about, we first cite the classical 
paper [17] and then [3, 6, 7, 19, 20] which treat problems close to our, without considering the two-scale 
aspect. 

Next we start to detail the method, by setting a generic and formal problem 

E'u'' = 0, (1.1) 
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in which E"^ is a partial differential operator depending on a parameter e which induces oscillations of 
typical size e in the solution u'^ = u'^{z). In the case where e is uniformly small over the domain on which 
(1.1) is posed, a Two-Scale approach may give good results. To summarize, this approach consists in 
noticing that, in some sense, 



(1.2) 



where U = U{z, () is the Two-Scale hmit of u'^ (Two-Scale convergence goes back to [23]), and that U is 
the solution of a well-posed problem of the form 



£:[/ = 0, 



(1.3) 



where £ does not induce high frequency oscillations. Building a Two-Scale Numerical Method consists 
then in solving, via an usual numerical method, model (1.3) and rebuilding u'^ using an algorithm based 
on (1.2). 

In the case where e is not a uniform parameter over the domain, in the sense that it may be small 
in some sub-region and of order 1 in other ones, Asymptotic-Preserving schemes are intended to be 
used, since they allow to simulate both sub-regions with the same method. The methodology for these 
numerical schemes is based on the fact that u'^ may converge, in a strong or a weak topology, towards a 
limit u which is solution of 

Eu^O. (1.4) 



Then, an Asymptotic-Preserving scheme consists in solving 
such that 



and 



'■Az 



'■Az 



UAz 



as Az 0, 



as e — 0, 



where uaz is a numerical approximation of u, meaning that 

UAz — u , as Az — 0, 

and is the solution of 



Eaz UAz 



0. 



(1.5) 

(1.6) 
(1.7) 

(1.8) 
(1.9) 



where E^, is a numerical operator. Notice that the property (1.6)-(1.7)-(1.8) is called the Asymptotic- 
Preserving Property. 

All what was just said may be summarized in the following diagram: 



u'^ solution of 



solution of 



E' 



Az ^Az 







u solution of 
Eu^Q 



Az^O 



MAz solution of 

Eaz UAz = 



(1.10) 



In the case when u'^ contains high frequency oscillations with high amplitude when e is small, solving 
(1.5) can be time consuming to capture the oscillations. In this case, the convergence 



u, 



(1.11) 
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occurs in a weak sense only. Besides this, the Two-Scale limit U = U{z, Q may describe very well the 
oscillations of in the sense that the convergence 



(1.12) 



occurs in a stronger sense than (1.11). This is particularly the case when E'^ generates pseudo-periodic 
oscillations, essentially at a sole period, of size e in the solution u'^ when e becomes small. Moreover, u 
and U arc linked by a relation of the kind 

uiz) = j^J^U{z,C)dC, (1.13) 

where Z is the cell such that C H> [/(z, C) is Z-pcriodic for any z, and \2\ is its measure. 

In this case, we propose to add a layer in the diagram (1.10) which summarizes the ideas on which 
Two-Scale Asymptotic-Preserving Scheme is based: 



m'^ solution of 
E^u^ =0 



e — > , weak-* 



A2->0 



, Two-Scale 




u solution of 

Eu^O 



U solution of 
£U = 




solution of 
5^ C/^ = 



solution of 



E% 



'■A: 



= 



Az-i-O 



uaz solution of 
Eaz uaz = 



Uaz solution of 

£az Uaz = 




U%^^ solution of 



^Az "^Az 







(1.14) 

Diagram (1.14) may be looked at as a prism. At the bottom of the front edge stands the TSAPS and 
the back rectangle is made of the Asymptotic-Preserving scheme diagram of (1.10). The intermediate 
layer gathers the Two-Scale limit problem, at the top, and the Two-Scale Numerical Method, at the 
bottom. 

The solution U^^ computed with the TSAPS must be close to m*^ when it is computed in C = ^ and 
when Az converges to 0. This is the first part of the Two-Scale Asymptotic-Preserving Property and it 
is symbolized by the arrow going from [/^^ to and by the one linking and w*^. 
The second part of the Two-Scale Asymptotic-Preserving Property is symbolized by the right part of 
the diagram (1.14) : [/^^ must be close to the numerical approximation Uaz of the Two-Scale limit of 

^ /■ [ NUM I 

Ue when e — s- and a numerical integration (symbolized by j— / of Uaz needs to recover an 

1-2^1 Jz 
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approximation uaz of the weak-* limit u of u'. 

We will deduce the Two-Scale Asymptotic-Preserving model on the base of a Two-Scale Macro-Micro 
decomposition that yields an equation which stands at the top of the front edge. Arrows linking U'^ and 
u*^, and U and u may be seen as the continuous Two-Scale Asymptotic-Preserving Property. They 
are constraints for the Macro-Micro decomposition. We refer to [3, 6, 18, 20] for the use of micro-macro 
decomposition in order to design an asymptotic-preserving scheme. 

The present paper implements the Two-Scale Macro-Micro decomposition in the case of the linear 
Vlasov equation with a strong magnetic field. This equation is the following: 



^ + V ■ V^r + (E' + V X (B' + _) ) . = , 
r(i = 0,x,v) = /o(x,v), 



(1.15) 



where = /'(t,x, v) is the distribution function of an ion gaz submitted to external electric and mag- 
netic field = E'(t,x) and = B%t,x), where t e R+, x e K^^ ^ g ]^3^ ^^^^j where X = 27rei, i.e. 
the first vector of the canonical base (ei, e2, 63) of M^. 

The paper is organized as follows: in section 2, we recall the results of Frenod & Sonnendriicker [12] 
giving the Two-Scale convergence of . Roughly speaking they claim that, when e is small, 

/"^(^jX, v) is close to G(i,x,u(-, v)) in a strong sense, (1-16) 
/'(t,x, v) is close to /(t,x, v) in a weak sense, (1-1''') 

and the equations for G and / are known (u(t, v) is defined by (2.9)). In particular the v-dependence 
of / is only through = v • ei, v± = a/ (v • 62)^ -I- (v • 63)^ and not through the angle a, which is such 
that v ■ 62 = v± cos a. Based on this fact, in section 3, we build a classical Macro-Micro decomposition 
of (1.15) by following the same ideas of Lemon & Mieussens [20]. It consists in decomposing in the 
following way 

/^(t, X, v) = m'^{t, X, v\\,v±) + n'^{t, x, W|| , a) (1-18) 

with € Ker((v x Ai) ■ Vv) and n'^ G Im((v x A4) ■ Vv). Indeed, we will see in details that this 
decomposition is unique so that it is valid for every e, small or not. Note that to"^ does not depend on 
a, as /, while n'^ does. Equations for m'^ and n'^ are set out, making up the Macro-Micro model. The 
drawback of this model is that, because of (1.17), the exact equality (1.18) implies that function n'^ is 
still highly oscillating, as To avoid this drawback, we develop the Two-Scale Macro-Micro framework. 
This is done in sections 4, 5, and 6, where we build the Two-Scale Macro-Micro approximation of f^. 
First, based on results of Frenod, Raviart & Sonnendriicker [10] we improve (1.16) by using a first order 
approximation, claiming that when e is small 

/'^(t,x, v) is close to G(i, x, u(-, v)) + e Gi (t, x, u(-, v)) -f e /(t, -, x, v)) in a strong sense, (1-19) 

where Gi is the solution of a partial differential equation and I is given by a formula (see [10]). Inspired 
by (1.19), we construct our Two-Scale Macro-Micro decomposition in the following way 

r(t,x,v) =G(t,x,u(J,v)) + eGl(t,x,Lj(J,v)) +eZ(t, J,x,v)+e/i^(i, J,x,v), (1.20) 

where G| is intended to be close to Gi when e is small and /i"^ is the corrector to be taken into account 
when the order of magnitude of e is 1. In other words, we use the Macro-Micro decomposition in order 
to make precise for every e the approximation in (1.19) which is valid only for small e. Then, we need to 
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show that the decomposition in (1.20) is unique. This is done by paying attention to the properties of 
the operator 

|-+(vx A^).Vv, (1.21) 
since the expression (1.20) may also be set in the following form: 

nt, X, v) - A'-{t, J, X, v) + B\t, X, v) , (1.22) 

with 

A' = A'{t,T,x,v) eKei-{—+{vxM)-V^) and = 6'=(t, r, x, v) e (Ker(— +(vxX)-Vv)) . (1.23) 

Then, we obtain the equations for the Macro part G\ and the Micro one /i*^. In section 7, we underline 
the continuous asymptotic-preserving property of the Two-Scale Macro-Micro decomposition. We close 
the paper with two appendices about some properties of the operator in (1.21) leading mainly to the 
uniqueness of the decomposition in (1.22). 



2 Two-Scale convergence of 

In this section, we briefly recall the results of Frenod & Sonnendriicker [12] where the asymptotic be- 
haviour of (1.15) is studied. What is known and taken for granted concerning (1.15) is that, considering 
the following hypothesis 

/o>0, foeL^iR^), (2.1) 

and 

E^^E strong in L-(0,T;LL(R3)), 
B strong in L°° (O, T; iL(K^)) > 

we claim that, for any T > 0, (/'^)e>o is bounded in L°°(0, T; L^(]R^)) . Moreover, as e — )• 0, up to a 
subsequence, (/*^)c>o Two-Scale converges towards the Two-Scale limit (see [2, 23]) 

F = F(t,T,x,v) eL-(0,r;i^^(M+;L2(R6)))^ (2.3) 

where the subscript 7^1 indicates that r i— > _F(t,T, x, v) is 1-pcriodic for any (t,x, v). Introduced in [23] 
and further analyzed and used in several homogenization problems in [2], the Two-Scale convergence 
means that for any regular function ip = i'itjTjXjv) such that (t,x, v) 1— > -0(t,T, x, v) is compactly 
supported in = [0, T) x for any r, and r 1-^ ^(t, r, x, v) is 1-periodic for any (i, x, v), we have 

lim / p {ipy dtd^(hf = I I FiPdTdtdxdv , (2.4) 



where 



(V')^ = (V')^(<,x,v)=V(t,-,x,v). (2.5) 

e 



It was shown in [12] that the Two-Scale limit F has the following property: 



d 

F e Ker(— + (v x X) • Vv) , (2.6) 
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where ^ + (v x A^) • Vv is seen as anti-symmetric, non-bounded, with open range from L°° (O, T; L^_^ (E+ ; 
L^(R^))) on itself. The direct consequence of (2.6) is that, introducing the rotation matrix r as 

/I \ 

r(r) = cos {2ttt) - sin (27rr) , (2.7) 
y sin {2t:t) cos (27rr) J 

the Two-Scale hmit expresses itself in terms of a function G = G(t, x, u) e L°° (O, T; L2(]R3)) by 

F(t,r,x,v) =G(i,x,u(T,v)) , (2.8) 

where 

u(r,v) = r(r)v. (2.9) 

Using oscillating test functions defined by (2.5) for regular functions ip which are moreover in Ker(^ + 
{•V X A4) ■ Vv); it may be brought out that G is the solution of 

dG 

— + U|| ■ VxG + (E|| u X B||) • VuG = , ^2.10) 
G{t = 0,x, u) = /o(x,u) , 

where the notation U|| — (u • ei) ei is considered. 

Having in mind the diagram (1.14), we can say at this point that equation (1.15) plays the role 
of E"^ u*^ = 0, and (2.8)-(2.10) plays the role oi £ U = 0. So, if we build a numerical method to approxi- 
mate (2.10), this numerical method, coupled with (2.8) would play the role of £az Uaz = 0. 

In [12], it was also proved that the weak-* limit / of is the solution to 

Of 



. V|| •Vx/+(E||+vxB||)-Vv/ = 

/(t = 0,x,v) = / /o(x,u(r,v))dT. 
Jo 



(2.11) 



This equation is gotten from (2.8)-(2.10) by performing an integration in r and by using the link between 
the weak-* limit / and the Two-Scale limit F given by 

/(t,x,v) = / i^(t,T,X,v)dT. (2.12) 

Then, equation (2.11) plays the role of = in the diagram (1.14) and a numerical approximation of 
it would play the role of E^z uaz = 0. 



3 Classical Macro-Micro decomposition 
3.1 Rewriting the weak-* limit model 

Before going further, we will rewrite model (2.11) in a form involving cylindrical variables in velocity. 
In this coordinate system, the weak-* limit is independent of the angle variable, and consequently it is 
easily gotten that the v x B|| piece in the force term in (2.11) is removed. 

Because of the average over t and of the definition of u(t, v), it is obvious to see that the initial data 
in (2.11) only depends on x, V|| and v± = \/ ^2 + ^^f ■ In other words, writing the velocity in the following 
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cylindrical variables 

U|| G R such that U|| = v • ei , 



v_i_ e ]R+ such that v±_ = + , (3.1) 
a e [0, 27r] such that V2 ~ u_l cos a and = v±^ sin a, 



there exists a function mp = too(x, such that 

mo(x, V|| , = — / /o(x, u(r, U|| ei + wj_ cos a 62 + I'd, sina 63)) da . 



(3.2) 



The quickest way to get this is to notice that 



"II 

u(t, Gi + fj_ cosae2 + sinaes) = I cos(q; + 27rT) | , (3.3) 

v± sin(Q; + 2ttt) 

and that, consequently, integrating in r suppresses the a-depcndcncy. Hence, (3.2) writes also 

1 f^^ 

mo{x,vii,v±) ~ — / fo{x,vii ei + v± cos a e2 + v± sin a 63) da . (3.4) 
27r Jo 

Concerning the function / itself, we introduce the function m = m(t, x, , u^, a) linked with / by 

m{t, X, U|| , f_L, a) = f(t, X, W|| Gi + vj_ cos a G2 + u± sina G3) . (3.5) 

Because of (2.12) and (2.8), and using the same argument as just above, it is easily seen that m does not 
depend on a and that 



m(i, X, , w_l) = / G(t, X, W|| Gi + cos(q! + 27rr) G2 + w_L sin(a + 27rr) G3) dr 
Jo 

1 f^'' 

= — / G{t, X, VI iei+v± cos a G2 + v± sin a 83) da . 
27r Jo 



(3.6) 



Noticing that 



dm , . df 

— — [t, X, v\\,v±) = cos a — — (t, X, vm Gi + vj_ cosa G2 + v± sina G3) 
ovj_ av2 

df , . , 

+ sina — — (t,x, wii Gi + v±_ cosaG2 + w_l sin 063) , 

du3 

dm s df 

— — (i, X, wii , vj_) = —vj_ sina — — (t, x, vuei + v± cos ae2 + v± sina G3) 
oa OV2 



(3.7) 



df ^ 

+ v±_ cos a — — (t, X, ui I Gi + v±_ cos ae2 + v±_ sin a G3 ) , 
0^3 ' ' 

the equation for m can be deduced from (2.11). For this purpose, we develop v x B|[ in the new variables: 
u_L(cosaG2 + sinaG3) x B|| = |i?|| | (sinaG2 — cosaG3) . (3.8) 
As a consequence, we have 

(t;_L(cosa G2 + sina 63) x B||) • Vv/(i, x, t;|| Gi + v± cosa 62 + v± sina G3) 

dm^ , (3.9) 

meaning that the force term in (2.11) may be reduced to E||. Then, m is the solution of 



dm dm dm 
dt " " 
m{t = 0,x,v\\,v±) = moi-K,v\\,v±) 



9^ + ^119^+^119^^'^' (3.10) 
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3.2 Macro-Micro decomposition 

Inspired by (2.11) or more precisely by (3.10), we now write a classical Macro-Micro decomposition of 
the solution of (1.15). The aim of this decomposition is (see [3, 20]) to lead to a new model which 
will be equivalent to (1.15) when e ~ 1 and will formally give the weak-* limit model (3.10) when e — >■ 0. 
It is evident that this last property is not satisfied by the model (1.15). The new model will be the base 
for the build of an Asymptotic-Preserving scheme which will have to satisfy diagram (1.10). 

The main idea of the Macro-Micro decomposition is to write as a sum of a function being in 
Ker((v X A4) ■ Vv) and a function being in (Kcr((v x A4) ■ Vv)) . Adapting this idea to our context and 
having in mind (3.10), we seek mf = mf (t, x, U|| , v^) e L°°(0, T; L^(M''' x M x M+;t;xdxrfw|| duj_)) and 
p' ^ p'{t,x,vn,v±,a) e L°°(0,T;L|,JM;L2(M3 x R x m.+ ;vj_dxdvndvj_))) such that 



(3.11) 



and 



f^it, X, U|| ei -t- v± cos ae2 + v± sin a 63) = TO(i, x, U|| , wj_) + ml(t, x, W|| , uj_) + n'(i, x, wn , a). (3.12) 

The subscript ^27r indicates that a 1-^ p'^(i, x, V|| , u^, a) is 27r-periodic for any {t,x.,v\\,v±). 
Considering the anti-symmetric, non-bounded operator with closed range^ 



d 



L°°(0,r;i|^JI 



:L^(R^ X 



x M^; dxdv\\ dv±))) 



(3.13) 



we clearly have that 



e 

G 

G 



Then, if we inject the decomposition (3.12) in Vlasov equation (1.15), we obtain 



dm 
'dt' 



dt 
f 




Vx(m 



j -I- v± (cos a i?| — sin a B^) 
cos a -I- sin a £'| — v\ | (cos aB^~ sin a _B|) 



— ( — sin aE2 + cos aE^ + v\\ (sina + cos 0-82)) ^ ^\\ 



1 



/ dm 



dv\\ dv 



^11 
dm 

dv± 



II 

dm\ 
dv± 

da 



dn" \ 

dv\\ 
dn" 

dv± 
) 



^See Appendix A 



= 0. 



(3.14) 



(3.15) 
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Combining this equation with (3.10), we finally obtain 



dml 









v± COS a 


• Vxm + 


v± cos a 


\ v± sin a J 




I v± sin a / 



1 

vj_ 



1 

v± 



E^^ — + v±_ (cosaSI — sinaSI) 
cos a i?| + sin aE^ — 7;| | (cos a — sin a -B|) 

sin a i?2 + cos aE^ + \ (sin a + cos a -B|)) - 



+ v± (cos aB^ — sin a 
cos a E'l + sin a — W| | (cos aB^ — sin a 

sin a i?| + cos a i?| + U| | (sin a i3| + cos a -Bl)) ~ -^f| 



/ dm \ 

V y 

9w 



dml 

dv I 



V 



da 



(3.16) 



By projecting (3.15) onto Ker(^) and (Kcr(^))''~ = Ini(^), we will obtain the Macro-Micro de- 
composition of f^. We refer to Appendix A and [5] to get that projecting a function onto Ker(^) consists 
in computing its average in a and projecting it onto Im(^) consists in substracting from it its average 
value. This is what we do in the next lines. 

Since m, m|, E, E'^, B and do not depend on a, and recalling the definition of n'^ given in (3.11), 
the projection of (3.15) onto Ker(^) gives 



dr 



dr 



dt 



dxii 



E\ 



dm\ 

dv\\ 



dm 



27r 






sin a — — da 

COS a 



1 

— E'= 
2tt 



2-n 



sm a y 
COS a 









dp" 

dv± 
dp^ 

^dv± 

2tx 







da 



dx2 
dp'^ 

dx3 



da 



Vl 



—) da 

Vl' 



(3.17) 




f , , dp' 

/ (cosQ!i?2 + sinaBo) — — da 

Jo _ 



V 



2-K 



COS a I 



'•dv\^ 
, dp" 
^dv±_ 



—) da 
— ) da 



J 



In the second, third and last terms of the right hand side of (3.17), we use integrations by parts to tackle 
a-integrations. Finally, we perform the projection onto Im(^) by substracting (3.17) from (3.15), we 
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obtain a second order equation for p'^ in which m and m\ appear as right hand side members: 



dt da 



V 



cos a 



d^p' 

dx2 da 
d'^p' 

dxR da 



d^p' 

dx\ \ da 
1 

2^ 
1 

2^ 



d^p' 
cos a — -— da 



sin a 



dx2 da 
d^p' 

dx3 da 

d^p' 



da 



cos a 



1 

dv± da 27r 
d'^p' 1 

dv±^ da 271 
( 



-v\\ Bl 



1 9V 



dv\ \ da 

cos a — -— da sin a „ 

dv± da v± da'^ 

sin a — — da H cos a - 

dv± da v± da'^ 

d^p' 

dv\ \ da 27r 



27r uj 
1 

27T V 



' . d'^p" 
sin a „ ^ da 



2ix 



COS a 



cos a 



sma 



(cos aB\— sin a B\ 
d^p" 



dv± da v± 
d^p' 



1 . d^p 
sin a 



1 



i9u_i_ da v±_ 



cos a 



d^p' 
da^ 



1 

2^ 



1 /•2- 

— / (cos aBl — sin a B. 
d'^p" 



do? 
d'^p" 

da^ 

d'^p' 



da 



da 



-/( 



cos a 



sma 



dv± da 
d^p' 

dv± da v± 



^ dv^da 

1 . aVx 

— sin a 

vj_ 

1 d''p- 
cos a 



da'^J 




• Vxm - 



\ 

dm 



cos a 



dm 



sma 



(cos a i?| — sin a -B|) 



9w± / 
/ dm 
dvii 



vii B% 



cos a 



dm 

dv± 
dm 




(cos aB^ — sin a 

-fii SI 







\ 

dml 

dv± 
dml 

dv± / 

/ dm\ \ 



55 



dm\ 

dv± 
dm\ 

dv±_ I 



da 
da 



(3.18) 



10 



Concerning the initial conditions for m\ and p*^, we invoke the initial distribution /o as follows: 
/o(x, U|| ei + v± cosae2 + v± sin 0:63) 

1 r^- 

= — / /o(x, wii ei + v± cos + v± sm9e^) dO 
+ /o(x, W|| ei + uj_ cosae2 + w± sinaes) 



— — / /o(x, wii ei + v± cos 6* 62 + sin 6* 63) dO 
+ [/o(x,W|| ei + uj_ cosae2 + v±_ sinQ;e3) - mo(x, W|| , 



Since 



mo, we have 



9pe 



mi(i = 0,x,V||,u^) + — = 0,x,?;||,?;j^,a) 
oa 



= /o(x, W|| Gi + cos a 62 + w_L sin a 63) — mo(x, V|| , . 
By integrating this equation in a, we find 

■m\{t = 0, X, W|| , w^) = , 

so we also have 



p'^{t = 0, X, W|| , w^, a) = / /o(x, f II ei + v±^ cos 6 e2 + v±^ sin 63) d6 ^ a mo(x, v\\,v±_) . 



(3.19) 



(3.20) 



(3.21) 



(3.22) 



Remark 3.1 (Conceiving Asymptotic-Preserving schemes). Building a somehow discretization of 
the coupled problem (3.10), (3.17), (3.21), (3.18), and (3.22) would give a scheme which will be consistent 
with (1-15) when e > and consistent with its weak-* limit (2.11) when e is small. This scheme would 
play the role of the problem E^^u'^^^ — in the diagrams shown in (1.10) and (1.14). 



4 Two-Scale Macro-Micro decomposition: preliminaries 

The aim exposed in the first paragraph of subsection 3.2 will now be changed in order to obtain the 
Two-Scale limit model (2.10) when e ^ 0. Eventually, we will see in section 7.2, that the Two-Scale 
Macro- Micro decomposition described in the following sections, will lead by integration in r to the classical 
Macro- Micro decomposition. 

We assume fi'om now that the input electric field and the input magnetic field do not depend on e, 
i.e. we take 

E-^^E, B'=B. (4.1) 



4.1 Complements on convergence results 

Before entering in the core of the Two-Scale Macro-Micro decomposition, we integrate some knowledge 
which can be deduced from Frcnod, Raviart & Sonnendriickcr [10] where the asymptotic expansion of 
is presented. From the hypotheses (2.1) and (2.2) and adding some regularity assumptions for E and B, 
we can claim that 

1 / * \ 

- (^/^(t, X, v) - F(t, -, X, v)j Two-Scale converges to i^(i,T, x, v) as e — 0, (4.2) 
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and that F expresses in terms of other functions G and I in the foUowing way: 

F{t, T, X, v) = G{t, X, u(t, v)) + l{t, r, x, v) , 

where 

lit, T, X, v) = (r(T + - r(i)) V • Vx,G(t, x, u(t, v)) 

<r + ])- r{])) E{t, x) + u{t, v) x ((r(T + i) - r(i)) B(t, x) 



Vu^G(t,x,u(r,v)) 



and where G satisfies an equation of the form 



(4.3) 



(4.4) 



dG dG ~ 

■ — + (E,, + u • B||) • VvG = RHS{t, x, v, E, B, G) , 



at 



II 



(4.5) 



which right hand side can be exphcitly computed hke that in Theorem 4.2 of [10]. 

In these equations, G is the solution of (2.10), r is defined in (2.7), and the gradients Vx^ and Vu^ are 
defined as 

/ \ / \ 



Vx, = 



d 

8X2 

d 



Vu, 



d 
d 



(4.6) 



4.2 Sought shape of 

As announced in the introduction, we want to use a shape given by (1.22)-(1.23) to 
We clearly have that operator 



d 



+ (vxA^).Vv : L^{0,T;L^^iR+;L^{R^))) ^ L°° {0,T; L^^{R+; L^{R^))) , (4.7) 



is anti-symmetric, non bounded and with closed range (see Appendix B) and satisfies 

d d 
Kcr(— + (v X A^) . Vv) ® Im(— + (v x M) ■ Vv) = L'^(0,T;L^(R+;L^iR'^))) 

On another hand, we have 



(4.8) 



i^eKer(— +(vx7W)-Vv). 
Furthermore, wc can easily prove that the function I defined in (4.4) is such that 

d 

I e Im(— + (v X M) ■ Vv). 



(4.9) 



(4.10) 



Indeed, any regular function in Kcr(^ — h(vxA4)-Vv) reads 7(i, x, u(t, v)) forsome7 G _L°° (O, T; i-^(M^)) . 
Multiplying /(t,T, x, v) by 7(t, x, u(t, v)), integrating over x [0, 1] and performing the change of vari- 



ables 



V i-^. u = u(r, v) 



(4.11) 



12 



we obtain 

/ 7(tx,u) 



r{r'^)~r{-)dTj uj •Vx,G(<,x,u) 
^\-(r-i)-r(i)dr^ E(f,x) 



(4.12) 



r(T-i)-r(i) dr] B(t,x) 



Vu^ G{t, X, u) > dx du = , 



because the integrals over [0, 1] all vanish. 



Now, we integrate every remark we just did. We also notice that a decomposition of the type (1.22) 
cannot be unique because of choices concerning the place where oscillations are put: in the superscript e 
t . . 

or in variable -. This will also prescribe what we set now: we look for 



Gl = Gf(t,x,u) e L°°(0,r;L2(M6)) ^ 
= k'{t,T,yi,v) e L°°(0,r;L|'^(R+;L^ 

such that, for any r, 

fit, X, v) - G{t, X, lj(t, v)) + e Gl {t, x, u(t, v)) 

(dk^ \ 
—{t, r, X, v) + {vxM)- V^k^{t, r, x, v) 

From now, we shall denote 

(GoLj)(t,r,x,v) =G(t,x,Lj(r,v)) , 
leading, for instance, to the writing of = G o u. With this notation, we can claim that 

d 

(Gou) + e(Gl ou) e Ker(— +(vxA^)-Vv), 
e/ + e(^— +(vx7W)-Vvfc^) e Im(— + (v x 7W) ■ Vv) , 



(4.13) 



(4.14) 



(4.15) 



(4.16) 



and, as a consequence of (4.8), decomposition (4.14) exists and is unique since G and I are already iden- 
tified as the unique couple satisfying (2.10) and (4.4). 

In the previous section, when building the usual Macro-Micro decomposition, we projected the equa- 
tion satisfied by onto the kernel and the range of the involved operator. Here, this cannot be done. 

Indeed, here, the involved operator is — — I- (v x A^) • Vv and the function which has to be used in equation 

OT 

(1.15) in order to extract some informations on Gf and /c'^ is 

fit, X, v) = (G o u)(t, J, X, v) + e {Gl o u)(t, ^, x, v) 

t f dk^ t t \ 

+ e l{t, -, X, v) -f- e I — (t, -, X, v) + (v X X) • \7^k'{t, -, x, v) 1 , 



(4.17) 



where r was replaced by - in (4.14). Once this replacement done, we face with functional spaces where 

e 

operator (4.7) makes no sense. In order to bypass this difficulty, we use a weak formulation of (1.15) with 
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oscillating test functions 



i4>Y = (V')'(i,x,v) = V(i,-,x,v), 



(4.18) 



where ip = ip{t,T,x,v) is regular, 1-periodic in r and with compact support on [0,T) x M'' x M'' for any 
fixed T e [0, 1]. Writing 



(F^r = (Fi^)^(t,x,v) =Fi^(t,^,x,v) = (Gl ou)(t,^,x,v), 



(4.19) 



and using the same convention for other functions, the weak formulation with oscillating test functions 
reads 



\Fy+e{F^r+e{ir+e{hr]\{^y + ^{^y+^-{^.^r 



M. 



(E + V X (B + — )) • (Vv^) 



e ^ dr' 



dx rfv dt 



(4.20) 



/o(x, v) -0(0, 0, X, v) dx dv , 



where F is hnked with G by (2.8), F^ with by (4.19), I is defined in (4.4) and h" is linked with k' by 



OT 



(4.21) 



The way to deduce the equation for G\ and for k'^ consists in using, in a first computation, test functions 
lb G Ker( — + (v X M) ■ Vv), i.e. of the form 



(4.22) 

3 ^ Tn>3 Ti^Qj^^ a second compu- 



for regular functions 7 = 7(i,x, u) compactly supported in [0,r) x K'' x 

d 

tation, test functions belonging to — ^ x M) ■ Vv), or writing 

V' = ?^+(vxX).VvK, (4.23) 

OT 

for some regular functions k = /t(t,r, x, v) with compact support in [0,r) x R'^ x E'^ for every r S [0, 1], 
will be chosen. 



5 The Two-Scale Macro equation 



Wc begin by going further in the writing of weak formulation (4.20) by using the link between G and /. 
We recall this weak formula: 



Jvs- 



+ e 



+ e 



+ e 



{FY 



JR6 



+ V Vx(V')' + (E + v X (B + — )) • Vv(V')'J dy^dwdt 

, inr i^r + -i^r + v • (v^^)^ + (e + v x (b + — )) . (v^vo^ 



(0^ 



at 

dm 

dt 



dx dv dt 



M 

+ (vx — )-Vv(^)^ 



dx dv dt 



(ly V • (VxV^)' + (E + V X B) ■ (VvV^)* 



dt dx dv 



Jm<^ 



[ dih I dip , M , 

+ 7 + ■ + (E + V X (B + — )) . (Vv^)^ 



cZx dv dt 



/ /o(x,v) V(0, 0,x,v)dxdv. 

JK6 



(5.1) 
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Performing an integration by parts on the first term of (5.1) gives 

'-'^ <■ vdiibY , 7W , 1 

(^)' 4s^+v-Vx(V')' + (E + vx(B + — ))-Vv(V')' dxdvdt 
t L ot e 

^ I" \ dF \ dF M 1 

i^Y + 7 (7^)' + V • (VxF)^ + (E + V X (B + — )) • {V^FY (V')^ dxdv 



dr' 



/ F(0,0,x,v)?A(0,0,x,v)dxdv 

JR6 



(5.2) 



"'K" 



V • ( Vxi^)' + (E + V X B) • ( Vvi^)' (V')' dx dv fit 



/ /o(x,v)?/'(0,0,x,v)dxdv. 

JR6 



since _F G Ker(— — h (v x A^) ■ Vv). We invoke the hnk between F and G, and we consider the change of 

OT 

variables 



V n> u = u(-, v) which is equivalent to v = u( — , u) = u ^(-, u). 
e e e 

We also use the following notation 

(■0 ou"^)(i,T,x,u) = ■0(^,T,x,u~^(r,u)) = 7/'(i,T,x,u(-T,u)) . 
Since a direct computation yields 



(5.3) 



(5.4) 



iV^FY = V^{FY = Vv(G o u)^ = r(--)((VuG) ou) , 



(5.2) reads 



{FY 



^ + V • Vx(V')^ + (E + V x (B + — )) . Vv(0)^ 



dx dv di 



^ + {ri~i) u) . VxG + (r(^) [E + (r(-^) u) x B] • V.G 

(■(/) o I 

/o(x, v) ipiO, 0, X, v) dx dv 



' dt dx du 



(5.5) 



(5.6) 



^+ (r(--)u) • VxG+ [r(-)E + ux (r(-)B)] • VuG] (^ ou-^-^ dtdxdu 
ot e e e J 



- / /o(x,v)?/>(0,0,x,v)dxdv. 
Using now equation (2.10) satisfied by G, we finally obtain 



[FY 



^ + V • V^ii'Y + E + V X (B + — ) • Vv(V)^ 
Ot e 



dx dv di 



(r(--) u) . Vx,G + [(r(-) E)^ + [u x (r(-) B)] ^ 



Vu^G 



(5.7) 



/ /o(x,v) V'(0,0,x,v)dxdv. 

JR6 



(-0OU ^Ydtdxdu 
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Notice that the last term containing /g will vanish with the right hand side member of (5.1). 
Now, we study the third term of (5.1). For this purpose, we first notice that 

^i^^ ={i^)ou-^) + ((r(-r) u) X X) . ((VvO o u'^) 

= -(r(-T)u) -Vx^G - [r(T)E + ux (r(T)B)] -Vu^G, 



(5.8) 



hence, integrating by parts the third term of (5.1) and making the change of variables defined by (5.3), 
we get 



cZx dv dt 
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e{^y + {^Y + ivxM)-iS7jy {^Ydicd^dt- / (0'(0,x,v)(^)^(0,x,v)dxdv, 
e ( ^^^^ -) {^ou-^Y dx du dt 



JR6 



'ou ^Ydxdudt 



lo Jm8 



dt 

(r(--)u) • Vx^G+ [r(-)E + ux (r(-)B)] -S/^^g] Oj ou-^ djidudt 



(5.9) 



since Z(0,0,x, v) — (see (4.4)). The last term of (5.9) will vanish with the first term of the right hand 
side member of (5.7). 



Hence (5.1) has to be replaced by 
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dt 



[ip ou Y d^i du dt 



V dib 1 dih , M , 

■ + 7 + ^ ■ ^^'^'^^^ + (E + V X (B + — )) . (VvV)^ 



(0' v.(Vx7^)^ + (E + vxB).(Vvi^)^ 



dx dv 



(5.10) 



dx dv di 



{h'Y [{^Y + - (^)^ + V • (VxV')^ + (E + V X (B + — )) • (Vv^)^ 



dx dv dt = . 



Now, choosing ?/! as in (4.22) yields 

{Gl ouY ((-f) o u)' + V • ((Vx7) o u)' + (E + V X B) • Vv(7 o ■ 
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dt' 
dt ' 



dx dv di 



"f dxdudt 







^ 5(7 ou) 

^ at > 



dx dv dt 



(5.11) 



JREi 



(e (0^ + e {hy) V • (Vx(7 o u))' + (E + v x B) • (Vv(7 ° u))' 



dxdv dt ~ . 



Making in the first, third and fourth terms the change of variables v ^ u(-, u), replacing h'^ by its 

e 
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expression in terms of k'^ and dividing (5.11) by e finally gives 



^1 7? + (^(— ) ") • + E + u X (r(-) B)] . Vu7 



dx du dt 







9t 



7 dx du dt 



"'K« 



(Q°0^ + (K-5)-M)-((Vvfc^)ou-i)' 
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dx du dt 



JR" 



(lou-'r + {{—)ou-'y + {r{—)M) ■ {{w)ou-y 



X 



dr 

r(--)u) • Vx7+ [K-)E + ux (r(-)B)] ■ Vu7 



e ' " e 

which is the Two-Scale Macro equation of the model. 



dx du dt = , 



(5.12) 



6 The Two-Scale Micro equation 

In this section, we use in the weak formulation (4.20) oscillating test functions defined by (4.23). The 
computation leading the formula (5.10) is valid for any oscillating function -0 so we use it as a starting 
point for finding the Micro equation. 

d 

Recalling that Fi £ Ker(— + (v x M) ■ Vv), the second term in (5.10) yields 



Jm.<^ 



{Ftr [{^Y + -{^y + V • (V,V)^ + (E + V X (B + — )) • (VvV)^ 



dx dv dt 



^ r [ [Fir 

Jo Jr^ 



^+(vx— ).Vv(V.)^ 



dx dv dt 



V • (VxV^)'= + (E + v X B) • {V^ipy 



dx dv dt 



— e 



{—^Y{ijYd:>idvdt + e / ^^'(0, 0, x, v) ^(0, 0, x, v) dxdv 

;6 at 7^6 



el I {F^Y V • (VxV')' + (E + v X B) • (Vv0)' d^dvdt 
lo Jr6 L J 

ndF^ r Ok 
i^Y ij-Y + ^ M) ■ [v^^^Y d^d^dt 
,6 dt I dr J 

+ e / Fi'(0,0,x,v) — (0,0,x,v) + (v X yW) • Vvk(0,0,x,v) 



dxdv 



^ I I (FD' 

lo JR6 



v.(^)^ + (vx>().((VxVv«;)^v) 

+ (E + v X B) • i^^Y + {ExM+vxBxM)- {V^kY 

OT 



(E + v X A^) • ((V»^(v X A^)) 



dx dv dt . 



(6.1) 
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The first term of (5.10) reads 



e (— -) {ip o u-^Y dx dv dt 



dt 



fl 




JO JM 


6^ dt 'I 



dn t 
(— ou-i)^ + [(r(--)v) X M] ■ ((Vvk) ou-i) 



dx dv dt 



Tlic third term gives 



(ly V • (VxT^)^ + (E + V X B) • (VvV^)^ 



dx dv dt 



(0' V • + {vxM)- ((VxVv«)^ v) + (E + V X B) • 







+ (ExAl + vxBxAl)- (Vvk)" + (E + V X B) ■ ((V^^k)' (v x M)) 

Concerning the last term of (5.10), we have 
, dib 1 dtp , 

ihrii^r + -{^r)d^d.dt 

{h^Y —^d^dwdt 



dx dv dt . 



dt 



I I ---—{tpYdyid\rdt+ [ /i'=(0,0,x,v)^/'(0,0,x,v)dxdv 
Jo Jmb dt Jjj6 

'■'^ '■ 1 dh'^ dn 

+ / /i'(0,0,x,v) [—(0,0, x,v) + (v x 7W) • Vvk(0,0,x,v) dxdv 

i^Y + (V X A.) . (^)^ + -i^Y + -(V X A.) . (^)^) 



/O JRfi 



3k 

X [(^)' + (v X 7W) • (Vvk)'] dxdvdi 
or 



+ / [— (0,0,x,v) + (vx7U)-Vvfc'(0,0,x,v)l 
Jb6 ot 



rdn 



[— (0, 0, X, v) + (v X X) • Vvk(0, 0, X, v)] dx dv , 



on the one hand, and 



{h^Y V • (Vx^)^ + (E + V X (B + — )) • (Vv^)^ 



dx dv dt 



((^)' + (v X Af) • (Vvfc')^) [v- (^)^ + (V X A1) . ((VxVvAt)^v) 

+ (E + V X B) • (^^Y + (E X X + (v X B) X X) • (Vv^)'^ 
or 



1 



+ (E + V X B) • ((V»^ (v X X)) + - (v X X) • (— ^)^ 
+ ((v xM)y.M)- {V^nY + (v X M) • ((V^2ac)'(v x M)) 



dx dv di , 
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on the other hand. Remarking that 

((v xM)xM)- {V^kY = -(v X M) ■ ((Vvk)' x M) = ■ {V^kY 
and dividing by e, we finaUy get the weak formulation of the Two-Scale Micro equation: 

[(1^)' + (v X A^) • (Vvk)'] dxdvdt 
[—(0, 0, X, v) + (v X • Vvfc^(0, 0, X, v)] 



(6.6) 



Jmf' 



[—(0, 0, X, v) + (v X X) • Vvk(0, 0, X, v)] dx dv 
ii-g^T + (v X X) . (Vvfc^)^) V • i^Y + (v X M) . ((VxVvA^)^ v) 



dr 

(E + V X B) • (— + {ExM + {^rxB)xM)■ {W^kY 

OT 



1 



(E + V X B) • ((V»^ (v X M)) + -(v X A^) . i-^) 



Jki3 



(Vvk)^ + (v X X) • ((V»'(v X 

a7 



dx dv dt 



dx dv (ii 



(6.7) 



r 9V K (9V K 

[lY V • (^)^ + (v X A^) . ((VxVvk)^ v) + (E + V x B) • (— ^)^ 

+ (E X X + (v X B) X X) • (Vvk)' + (E + V X B) • ((V^^k)*^ (v x X)) 

dx dv dt 



dx dv dt 



(-^r (i^r + (V X A^) • (v.K) 







at 



F^iO, 0, X, v) —(0, 0, x, v) + (v X X) • Vvk(0, 0, x, v) 



dxdw 







V ■ {^Y + (v X A/() • ((VxVvk)^ v) 



(E + V X B) {^^^Y + (E X X + (v X B) X X) • (Vvk)"^ 
or 



(E + V X A^) • ((V»^(v X A1)) 



dx dv di = . 



Remark 6.1. We are aware that the Two-Scale Macro-Micro system, (5.12)-(6.7), has not necessarily 
a unique solution. In future work we will look for an additional condition for G\ and fc' leading to 
the uniqueness of the whole solution (Gf,fc'). Such a condition will be usefull for the conception of the 
numerical scheme for the Two-Scale Macro-Micro model in order to be sure that when e ^ 1 we will 
approximate the solution of the starting Vlasov equation (1.15). 
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7 Asymptotic-Preserving Properties 



7.1 The Two-Scale Macro-Micro problem 

The Two-Scale Macro-Micro problem, i.e. the system of equations (2.10)-(4.4)-(5.12)-(6.7), is more 
complicate to solve than the original problem (1.15), but has two advantages over the latter. First, 
when e — !■ the Two-Scale Macro-Micro model reduces to the limit model in (2.10) while equation 
(1.15) becomes singular at the limit. Second, the Two-Scale Macro-Micro model decomposes the original 
solution into a macro part, G o u -I- e o u, and a micro part, el + eh^. This question is relevant 
when doing the transition from the very small e regime to the one of e ^ 1, since our model describes 
separately the evolution at the macroscopic time (of Gou and G\ ou) which contains essential oscillation 
through u, and the evolution of oscillation corrections [l and h^) at the microscopic time. 

Now it is easy to see that, under the hypothesis of uniqueness of the solution of equations (5.12) and 
(6.7) (see Remark 6.1), the Two-Scale Macro-Micro problem (2.10)-(4.4)-(5.12)-(6.7) is equivalent to the 
original problem (1.15). Indeed, let be the solution of (1.15). We have seen, by using Lemma B.l, that 
the decomposition (4.14) exists and is unique. Then, by the calculations we did in the previous sections, 
we obtain that G is solution to (2.10), I is given by the formula (4.4), G\ is solution to (5.12) and k"^ to 
(6.7). In particular, we have proved that the equation system (5.12)-(6.7) has solution. 

Conversely, assume that the solution of the system (5.12)-(6.7) is unique. Let (G, /,Gf,fc'^) be the 
solution of (2.10)-(4.4)-(5.12)-(6.7). Then constructing by (4.17), we wiU obtain the solution of (1.15), 
since this problem has unique solution. 

7.2 Macro-Micro Decomposition vs Two-Scale Macro-Micro Decomposition 

Property 7.1. Integrating in t the Two-Scale Macro-Micro decomposition leads to the classical Macro- 
Micro model developed in Section 3.2. 

Indeed, let us recall that the unique decomposition made in (4.14) is of the kind 

r (t, X, v) - F{t, T, X, v) + Flit, T, X, v) + Hlit, T, X, v), (7.1) 

for some Ff e Ker(^— + (v x M) ■ Vv) and HI G Im(^— + (v x M) ■ Vv) . Then, integrating (7.1) in t, 
we obtain using (2.12) 

r(t,x,v) = /(t,a-,v) + / Fl{t,T,x,v)dr + [ Ht{t,T,x,v) dr. (7.2) 
Jo Jo 

Now, changing the variable v in {v\\,v±,a) defined in (3.1), and using Lemma B.2, we obtain that the 
function (v^ , a) I— > / _F'2'^(<, r, x, U|| ei + w^(cos a e2 + sina 63)) dr is in Kcr^^— ^ . Similarly, {v±,a)i-^ 

/ Hf (t, T, X, U| I ei + v±_ (cos a 69 + sin a 63)) dr is in Im I — — ] . Next, writing (7.2) in the new variables 
Jo ^oaJ 

and recalling the decomposition of /"^(tjX, v) in (3.12), we deduce from (3.5) and the uniqueness of such 
a decomposition (stated in Lemma A.l(iii)) that 



TOj(t,x, W||,-y_L) = / F;[^(t, T, X, W|| Gi -I- t;_L(cos a 62 -I- sina 83)) dr, 
Jo 

n'^(i, X, W|| , u_L, a) = / iJj^(t, r, x, U|| ei + wj_(cos a 82 -I- sina 63)) dr. 
Jo 



(7.3) 
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7.3 Convergence of the Two-Scale Macro-Micro problem 

Let us recall that our Two-Scale Macro-Micro decomposition (4.14) is based on the convergence result 
in (4.2). More precisely, keeping in mind that G + eG + el is the first order approximation of f^, then 
(4.14) can be seen as a Macro-Micro decomposition at the first order level of approximation of f^. Next, 
we justify this approximation using the Two-Scale convergence. 

Theorem 7.2. We assume that /o G L^{R^), E G W^-°°{R^), B e W^-°°{R^), and dl/dt, Vx,v^ G 
L°° {0,T; L°°{0,1; L'^{R^))) . Then, when e ^ 0, the solutions G\ o u of (5.12) Two-Scale converges to 
G G L°°(0,r;i2(R6)^_ ^^jg solution of (4.5). When e 0, the solutions k" of (6.7) Two-Scale converges 
to 0. 

Proof Recall from Theorem 1.5 of [10] that when e -> 0, (^i (/' - G o u) - Two-Scale converges to 

G o u G Kcr(^— + {-v X M) ■ Vvj . Then, since (4.17) implies that 

^{f ~ G on) - I = Gl on + h\ 

and since G| o u G Kct(^— — h (v x M) ■ Vv^ and h'^ G — I" ('^ ^ ' , the theorem's conclusion 
is clearly true. □ 



Appendices 

In the first Appendix we characterize the projection onto Ker(^). We establish that it consists in 
computing its average in a and projecting it onto Im(^) consists in substracting from it its average 
value. Then, in the second one, we do the same with operator ^ + (v x A^) • Vv. We also establish the 
link that exists between these two operators. 



A About the operator — — 

da 



Lemma A.l. Let the unbounded operator —— : L\ (M; L^(R+; di;^)) L\ (K; L^(R+; dw^)). 

da """^^ 

Then, it has a closed image in L\ (M; L^(IR"'"; duj^)) and 



r^Ker(-) =lm(-). 



(a) For any function f G Li^^^iR; L? {W^ ; v ± dv ±)) , the projection of f on Ker( — j is the function 



Pf : {v^,a) ^ ^ r f{v±,e)dd. (A.l) 



(ill) L|^JR;L2(M+;„^rf„^)) ^Ker(— ) ©Im(^). 
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Proof, (i) Let q G ^Ker(^— )^ . Then, for any g e Ker^— ^, i.e. g = g{v±), we have 



27r /' + CXD 

q{v±, a) g{v±) v± dv± da = . (A. 2) 



/o ./o 

Remarking that q can be written as 

q{v^,a)^^(^j^ q{v^,e)de^, (A.3) 
we only have to prove that the function K defined by 

nC\. 

K{v^,a) = / q{v±,9)d9, (A.4) 

is in (R; i^(R+; dux)), i.e. is 27r-periodic in a and / / |/^(wx,a)| vj^dv^da < +00. 
First, we have 

r2ir 



K{v±,a + 27r) - K{v^, a) = q{v^,e)de = 0, (A.5) 

Jo 

by (A. 2). Secondly, using the Jensen inequahty, we remark that 

pa /'27T 

\K{v±,a)\^ <\a\ \q{v±,e)\^ d0 <2n \q{v_L,e)\'^ dO . (A.6) 



Then 

/ / \K{vj_,a)\^ v±dvj_da < Att^ / / \q{v±,9)\^ v± dvj_ d9 < +00 , (A. 7) 

Jo "'0 "'0 "^O 

since q € L^^^(K; L^(R+; dv±^)). Thus (7 e Im^— ^ and so Ker^— ^ C Im^— ^ . Since the converse 

inclusion is obvious, {i) is proved and the image of the operator is indeed closed. 

{n) It is clear that / G L'^^^^{R; L'^ {m+ ; v ± dv i^)) implies that PJ £ L'^^^^{R; L'^ {E+ ; v ± dv . Then, 
using the projection characterization, we have for any function ip e (R; i^(M+; wj^ dw^)) such that 

da ' 

^0 



f{v^,a)-— I f{v^,e)d9 ] ij{vj_)v^dvj_da 



+ 00 / /•2tt \ / i'2iT \ I r-2-K \ 



\ "'0 / \ Jo / Jo \ "'0 



= 0. 

/ d - 

Hence P/ defined in (A.l) is the projection of / on Kcr(^^— 

(ill) Using item (i) wc have only to prove that any / G (M; L'^(R+; duj^)) writes as a sum 

f = g + h with .g € Ker(^— ^ and ft, 6 Im(^— ^ Given / e L|^^ (R; gf^,^))^ let g be the 

projection of / on KerfT— ^ and let h ~ f — g. Then, by item (ii) we only need to show that h E Imf — — 

This can be done using a similar way as in the proof of item (i). This concludes the proof of the 
Lemma. □ 
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d 

B The operator — + (v x A^) • Vv 

OT 

Lemma B.l. Let the unbounded operator 
d 

■+(vxX)-Vv : L°°(0,T;L^^(M+;L2(]R6))) ^ L°° (O, T; (M+; L2(M6))) . Then we have 



dr 



L°^{0,T;L^^{R+-L^{M.'^))) = Ker ( + {^r x M) ■ e Im + (v x A^) • Vv j . (B.l) 



dr 

Proof. Let f e L°° (O, T; (M+; 2.2(^6))) . Through the change of variables 



U|| G K such that W|| = v • ei , 

v± e M+ such that v± = + w| , (B.2) 
a G [0, 27r] such that U2 = u± cos a and U3 = v±^ sin a, 



/ ^where /(f, T, X, U|| , a) = /(t, r, x, U|| ei + ^^(cosa 62 + sina 63)) j is in L°° (O, T; (I 

9 d d 

— + (v X A^) • Vv writes 27r— . 

OT OT Oa 



d d d 
L^fM^ X E X R+:v±dxdv\\ dv±)))) while the operator — + (v x M) ■ Vv writes 2tt—. Next, 



performing the second change of variables 



(7 = T+ —a, 

(B.3) 

B = T a , 

^ 27r ' 



/ (where/(i,a,x,t;||,«^,/3) - /(i, (a + /3)/2, x, w^, ^(a - /?))) stays in (O, T; (R+; ^^^^^(I 

d d , d 

— 2777— (and consequently — 

OT Oa OT 



3d 
L^(M^ x M X M+; w_L dxdv\\ dv±)))) and the operator 27r— (and consequently — + (v x M) ■ Vv 



d „ 
becomes From now on, let us fix (i,x) G [0,r] x R . Considering 

^ :i^^(R+;i|^^(R;L2(M3 xMx]R+;i.xdxdu||dw±))) 

i|'^(R+;i|^^(R;L^(R^ x R x R+;v^dxdv\\ dv^))), 

d\ / d 



we know, by Lemma A. 1, that L^^^{R;L^{R+;v_i_dv±_)) = Ker(^— j ©Im(^— j . Hence, because (R+ ; 

L|^JR;L2(k X R+;w_Lrff|| rff_L))) C (R+; JR; L2(R x R+; t;_L dw_L))), we can write / = 
Pf + if - P.f), where Pf is defined as in (A.l). Since / was chosen in L'^^{R+ ; L^^^JR; {M. x 
M.^;v± dv^^ dv±))) then it is easy to see that Pf belongs to this same space. Obviously / — Pf G 
L^^(R+;L|^^(R;i2(R x R+;vj_dv\\ dvj_))) and thus (B.l) is proved. □ 

Lemma B.2. Let f G L°°{0,T;L^^{m.+ ;L'^{M.^))) and let f be defined as in the proof of Lemma B.l. 
Then 

(i) if f £ Ker^^+ (v x M) ■ Vv^ then for any fixed (t,x, W||) G [0,r] x R^ x R we have 

T : (wj.,q:)i-^ y /(t,r, X, W||,w_L,a)rfr G Ker(^— ^ 
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(a) if f & lin\^— + (v X M) ■ Vv j then for any fixed {t, x,v\\) G [0, T] x x M we have 

"1 , a . 



T:{v^,a)^ I f{t, T, x,vii, v±, a)dT elTa(^—j. 



Proof {{) Let / e Ker(^— + (v x X) • Vv) and fix (t,x,W||) G [0,r] x x R. Thus, by the proof of 

~ / d d \ 

Lemma B.l, / e Kerl — 27r^:^) . We then derive J- 



V Ot daJ 
dT 



da 



(t)_L,a) 



19/ 

da 

1 r^df 



2tt Jn dr 


d 



(t, T, X, w_L, a) dr 
(t, r, X, wii , u^, a) c?T 



(B.5) 



by the 1-periodicity in r of /. Thus, F e Kerl -tt-)- {ii) Now let / e Im( (v x M) • Vv I . Then there 



da' 



d 



dr 



, , o ^ ^ ~ dh dh 

exists a function /i g i°° (O, T; (M+; £2(^6))^ such that / = - — Stt— . Fix (t,y.,v\\) e [0,T]: 

Integrating in t from to 1, we obtain, using tliat h is 1-periodic in r, 

d 



leading to £ Iml — — and thus concluding the Lemma. 
\da/ 



□ 
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